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The influence of jet geometry on light curves 
and spectra of GRB afterglows 
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Abstract. We have performed detailed calculations of spectra and light curves of GRB afterglows assuming that the observed 
GRBs can have a jet geometry. The calculations are based on an expanding relativistic shock GRB afterglow model where the 
afterglow is the result of synchrotron radiation of relativistic electrons with power-law energy distribution at the front of external 
shock being decelerated in a circumstellar medium. To determine the intensity on the radiation surface we solve numerically 
the full time-, angle-, and frequency-dependent special relativistic transfer equation in the comoving frame using the method of 
long characteristics. 
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At the present time we know that gamma-ray bursts (GRBs) are 
explosive phenomena at cosmological distances. If the emis- 
sion is isotropic, estimations based on observations give us 
the values of released energy up to Eq ~ 3.4 x 10 54 ergs for 
GRB990123, that exceeds the rest energy of a solar mass star 
(Kulkarni et al. 119991) . To reduce this large amount of energy 
it can be supposed that the GRB emission is highly collimated. 
May be, the better evidence for jet structure is the achromatic 
break in light cirves (Sari 1999) of the light curves seen in 
many afterglows, e.g. GRB990123 (Kulkarni et al. [19991 ) and 
GRB990510 (Harrison [1999l Stanek [1999l . And, finally, spher- 
ical symmetry conflicts with linear polarization (Sari 1999 ) ob- 
served for a few afterglows (Covino [l999l Wiiers [19991 ). 

Generally, a GRB jet can display an angular structure and 
can be seen by observer at wide range of viewing angles from 
the jet axis (Wei & Jin 2003, Granot 2003). For now, however, 
we consider a jet with uniform angular structure taking into ac- 
count the effect of equal-arrival-time surface at different angles 
of observation and show which changes in GRB afterglow are 
produced in transition from spherical symmetry to jet geome- 
try. 

The evolution of the jet and the light curves has been widely 
investigated (Panaitescu ll9991 Kumar 2000), including lateral 
jet expansion (Salmonson 2003), investigation of the 'stuc- 
tured' jet (Granot 2003 ), 3D numerical simulations of the jet 
dynamics (Canizzo 2003) and different angles of observation 
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(Granot 2002). All these works have not solve accurately trans- 
fer equation for the resulting light curves calculations and are 
based either on simple expressions for local emissivity or fo- 
cus on power-law branch of the spectrum between the break 
frequencies or some other simple assumptions on the charac- 
teristics of radiation field. 

In this paper we present a detailed calculation of spectra 
and light curves with preliminary numerical solution of special 
relativistic transfer equation in the comoving frame. We will 
show that the exact calculation of intensity, depending from 
the angle to the surface normal, can have an essential influence 
on the form of the spectra visible to the observer. The exact 
determination of equitemporal surfaces which is important for 
explanation of observed luminosity in GRBs and comprehen- 
sion of their spectral properties (Bianco & Ruffini [T976l) is also 
taken into account in our calculations. 

The transfer equation needs to be solved for the accurate 
calculation of intensity on the surface of radiating structure by 
integrating the emission along the characteristic through the 
structure and the following fluxes calculations. The calculation 
is based on the model where the jet is cut from a spherically 
symmetric flow. We add two parameters for taking into account 
a radial jet structure and different values of observer viewing 
angles. In the next section we discuss our model in more de- 
tail, in section 3 we calculate the emission for different values 
of the parameters and finally we present some discussions and 
conclusions. 
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Fig. 1. The part of quasi-ellipsoid from which the photons reach 
a remote observer. The explosion center is located at the vertex 
of the jet opening angle an m . The jet is observed under the angle 
a bs between the jet symmetry axis and the direction towards 
the observer. 

2. Physical Model 

Our jet geometry investigation is based on the numerical solu- 
tion of the problem in the spherically-symmetric case. To take 
into account jet geometry first we fix the direction towards the 
observer in the case of spherical symmetry. Because of the high 
shock velocity light at a certain time reaches the observer from 
the ellipsoidal structure, the part of which can be seen in Fig.Q] 
To construct the jet we cut the cone with axis being the di- 
rection towards the observer A and the angle a\i m , forming jet 
opening angle. The sight of observer A in this case coincides 
with jet axis. To consider the jet under different values of view- 
ing angle we add angle or b s between the observer B and the jet 
axis. This approach gives us the possibility to consider some 
jet effects without using specific hydrodynamics code and more 
complicated transfer equation. 

Let us consider the model we have used for numerical cal- 
culations of spectra and light curves in the case of spherical 
symmetry (Tolstov & Blinnikov 2003). 

In general, the transfer and hydrodynamic equations con- 
stitute a combined system of equations. In our problem, how- 
ever, we solve them separately. To determine the variables of 
the medium we use a self-similar solution for a relativistic 
shock for an ultrarelativistic gas in the case of the conserva- 
tion of total shell energy (Blandford & McKee 1976). The so- 
lution describes the explosion with a fixed amount of energy Eq 
and propagation of a relativistic shock through a uniform cold 
medium. 

1 ^ t^2 _17/12 2 l-r-2. -1 -> i-2 -7/4 

P=3 e =3 Wl T X ; T = 2 Tx ' n 7 = 2n i r X - 
where 

r 2 ocr 3 ; x = [i + 8r 2 ](i -r/t). 

Above, T is the Lorentz factor of the shock front, y, p,e,n - 
the Lorentz factor, pressure, energy and density of the shocked 
fluid, measured in the local rest frame of fluid, respectively, and 
«i is density of the external medium. 



For accurate calculation we should know also the electron 
energy spectrum and the magnetic field strength. Here for lo- 
cal emissivity calculation we use the conventional assumptions 
from relativistic electrons (e.g. Sari 1999 ). We assume being 
based on standard fireball shock model (Zhang 2003 ) that the 
electrons have a power-law distribution and that their total en- 
ergy behind the shock front accounts for e e of the internal en- 
ergy: 

N(y e ) = K y~ p ; y e > y min>0 = — , 

where m e is the electron rest mass and Kq = (p - Y)noy p ~ „. 

The magnetic field is parameterized by the quantity eg, 
which is equal to the fraction of the internal energy contained in 
the magnetic field B 2 = Snese. The magnetic field is randomly 
oriented and decreases with time due to adiabatic expansion of 
the shell. 

As the electrons pass through the shock they begin to lose 
the energy through adiabatic cooling determined by the solu- 
tion of Blandford & McKee ( I19761 ). This process is well de- 
scribed in detail by Granot & Sari (2002). 

Here we present only the basic formulas for synchrotron 
radiation used in our calculation. The spectral power of a single 
electron averaged over the pitch angle is: 

P{a>) = - -F(-), 

where 

1 i o 3n eB 1 

Psy = —cr Th cB 2 (y 2 - 1); oj c = — y 2 , 

on 8 m e c 

and F(u) - a standard synchrotron radiation function (Rybicki 
119791 ). The synchrotron absorption coefficient is specified by 
the formula: 

1 p"« N(y) d I 2 \ 
X %mn e v 2 J rmm 7 y 2 dyV ' / 

To determine the intensity on the radiation surface we solve 
numerically the full time-, angle-, and frequency-dependent 
relativistic transfer equation in the comoving frame using the 
method of long characteristics up to the values of Lorentz- 
factory ~ 1000. 

After the calculation of intensity the flux can be deter- 
mined: 

^0,^ = ^2- J /jR 2 I(r(p ), v (— J,cos(5(cos£o))^j dp Q , 

where subscript o is related to the observer frame, D is the dis- 
tance to observer, p = R/D. p mm - cosine of the maximum 
angle visible to the observer. 

To sum up, the observer afterglow spectra and light curves 
depend on the the hydrodynamic evolution, the radiation pro- 
cesses, the distance to the observer and the two parameters we 
have used for taking into account jet structure: the jet opening 
angle a\\ m and the viewing angle a b s . 

In our calculation for solving the problem in the case of 
spherical symmetry we have used the following parameters: 
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Fig. 2. The angle visible to observer from emitting structure 
increases with time. Time is measured in the observer frame of 
reference. For the t — 10 4 s it is shown that there is no reasons 
to increase the jet opening angle more that 0.037. 



Fig. 4. Afterglow light curves at frequency v = 5 x 10 14 Hz and 
observational angle a \, s = for different values of limitation 
angle ai im . 




12 14 16 18 

lg v, Hz, a = 



is the following relationship between the real observation angle 
9 and a has : 

1 - pp(p) 



[i+P 2 ( M )-2 P (p)p]^' 



where p. = cos a, p = R(t obs )/D, f obs = t obs (p, p, D), D - dis- 
tance from the center of the burst to the observer, and time in 
the burst frame of reference t is connected with the time in the 
observer frame of reference t \, s by the formula: 



tabs = t + 



D(\-[\+p 2 Qi)-2pUi)p\ l l 2 )+R 



Rq - R(tobs) , , \ 



Fig. 3. Instantaneous afterglow spectra at time t = 10 4 s and 
observational angle a Q b s = for different values of limitation 
angle a Um . 



E Q = 10 53 ergs, m = 1 cm" 3 , e e = 0.5, e B = 0.1, p = 2.5, 
D = 10 27 cm. In the next section we consider the results of the 
jet geometry influence on spectra and light curves varying a\\ m 
and a obs . 



3. Results of the numerical calculation 

As the shell from which the light reaches the observer moves 
towards the observer, the angle of the structure visible to the 
observer increases with time. This dependence is presented in 
Fig. |2] where a is the angle between the direction to the ob- 
server and the line connecting the center of the symmetry with 
the point on the surface that is still visible to the observer. There 



if we suppose that the initial time of observation corresponds to 
the initial time of the burst (to = tn, i,j + D/c) and initial radius 
of the burst structure is Rq. 

Now if we fix the time of observation by the value, say, of 
f', = 10 4 c, there is no reason to increase the value of lim- 

obs ' f 

itation angle in our jet structure for more than ar ]im = 0.037 
(Fig. 13, because this does not produce any effect on the result- 
ing spectra and light curves as if they are considered from the 
non-limited structure. 

In Fig. [3] and Fig. [4] we can see the calculated spectra and 
light curves at zero observational angle and at different values 
of limitation angle. 

The changed form of the spectra, having at some values 
of limitation angle two peak fluxes, is the consequence of 
the ring intensity structure on the radiative surface (Tolstov & 
Blinnikov l2003l . 

If we look at the ring structure we see that than more the 
light frequency than closer the maximum of brightness to the 
edge of the image. 
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Fig. 5. Instantaneous afterglow spectra at time t — 10 s and ob- 
servational angle a Q \, s = 0.01 for different values of limitation 
angle ai im . 



Fig. 6. Afterglow light curves at frequency v = 5 x 10 Hz and 
observational angle a b s = 0.01 for different values of limita- 
tion angle an m . 



The flux from the observed image can be calculated by the 
formula 

F v - 2n I /(cos 6, r, v) cos 8d cos 

Jcos 6 max 

where 9 is the angle between the point on the radiating surface 
and the center of the structure (see Fig. [TJ, 8 max corresponds 
to the edge of the image. In the absence of jet limitation angle 
the flux is monotonically decrease in increasing light frequency 
at the right part of the spectrum. In the presence of limitation 
angle some maximums of intensity at lower frequencies are ex- 
cluded from the flux integral and it gives at higher frequencies 
flux value compatible to that one at lower frequencies. 

The light curves do not show this effect just having "jet 
breaks" due to the limitation angle. This results from decreas- 
ing of the radiation arrived to the observer from the shock lim- 
ited by aiim. Larger the value of «/„„, less radiation at some 
frequencies gets towards the observer. Of order of days at 
au m = 0.2 we can see the break typical for some observed op- 
tical afterglows (Zhang & Meszaros l2003l l. 

In Fig. [5] and Fig. [6] we present the calculated spectra and 
light curves at the observational angle = 0.01. At small 
values of limitation angles we can see also the changes in spec- 
tra but for the light curves we did not have this effect. Some 
of the presented light curves are cut at early moments of time. 
This is the consequence of solving transfer equation only up to 
the Lorentz-factor value y = 1000. At larger value of Lorentz- 
factor we suppose that the matter is not radiative and if at 
ttobs = this effect is not revealed that at a D b s = 0.01 some 
of the characteristics start at Lorentz-factor value y > 1000. 

4. Conclusions 

It is widely believed that GRBs are born in jet geometry. In 
this case the resulting afterglow radiation becomes highly col- 
limated. The numerical calculations of light curves in these 



models (Granot 120031 Salmonson 120031 ) are based on some 
assumptions for intensity on the propagating shock front. The 
shape of the local spectral emissivity is approximated as a bro- 
ken power-law with some typical breaks corresponding to syn- 
chrotron radiation. As we can see from our results the spec- 
tra can have some peculiarities and the shapes different from 
power-law as in direct view to the jet as at some angle to the jet 
axis. 

Of course, our consideration does not take into account 
some effects of jet model and the exact numerical calculations 
should be at least two-dimensional to allow lateral expansion 
and angular structure of the jet. 

Nevertheless we would like to point out that the accurate 
calculation of intensity using relativistic transfer equation can 
have an influence on the shape of spectra and light curves of 
GRB afterglow. Constructing more precise model using the ex- 
act numerical calculations can help to explain the peculiarities 
of GRB afterglow and shed some light on the nature of GRB 
phenomenon. 
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